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ABSTRACT 

The destruction of planets by migration into the star will release significant amounts of energy 
and material, which will present opportunities to observational study planets in new ways. To 
observe planet destruction, it is important to understand the processes of how this energy and 
material is released as planets are destroyed. It is not known how fast the large amounts of 
energy and material are released, making it difficult to predict how observable planet destruction 
will be. There is a huge amount of energy made available by falling deep into the star's potential 
well: Simple calculations show that many of the currently known "hot Jupiters" can potentially 
produce events as luminous as a small nova if the energy is released fast enough. To observe these 
events, the important questions are how will this energy be released, and whether the energy will 
be released rapidly enough to create an event luminous enough to be found by transient surveys. 
Even the final rapid tidal infall of the planet may input enough energy into the star that could 
rival the energy output of the star. The final destruction will release even far more energy. These 
events are rare enough to explain why there is no undisputed observation of planet accretion, 
but future transient astronomy surveys such as PTF, Pan-STARRS, and LSST may have a good 
chance of catching planet destruction events. 

Alternatively, if planet destruction is slowed by the inward migration alternating with periods 
of outward migration caused by Roche lobe overflow (RLOF), then the primary signature may 
be the effects of the release of large amounts of gas. The infall of this gas also may significantly 
contribute to the system's luminosity. The release of planetary gas may be a searchable signature 
of planet destruction. Signs of runaway RLOF and outward or alternating RLOF should be 
searched for. RLOF unevenly distributed through the orbit may maintain an eccentricity through 
orbitally asymmetric transfer of angular momentum. 

Smaller planets may provide less dramatic but potentially more common events. Planet for- 
mation models frequently refer to protoplanet infall. Comparison with solar system objects shows 
how these events may be much brighter than seen in ordinary YSO and proto-star variability. 

Subject headings: Extrasolar planets, planetary dynamics, planetary evolution 

1. Introduction well understood, though much important work has 
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by observations of gas release through Roche lobe 
overflow (RLOF) that could be considered as the 
"beginning of the end" of the most close-in plan- 
ets (Fossati et al. 2010). The closest exoplan- 
ets now being found close to their stars may be 
tidally spiraling in to their destruction even dur- 
ing their main sequence. Though planet consump- 
tion has long been expected to happen late in the 
evolution due to explanion of the star, and has 
long been considered during planet formation, the 
earlier consensus was that such planets were not 
likely to migrate into the star during the main se- 
quence phase. Tidal migration has become one 
of the fastest developing areas of exoplanet study, 
with extraordinarily rapid exchange of new under- 
standings from new observational and theoretical 
results. Within the last few years there have been 
enough transiting exoplanets for patterns of tidal 
migration to become apparent through study of 
orbital parameters. New theoretical results and 
observations have brought new insights addressing 
the question of how the closest in planets get de- 
stroyed by tidally migrating into their stars. The- 
oretically, that the destruction of planets by in- 
spiral into the star actually must occur has been 
shown by recent work that considers the effect of 
tides induced on the star (eg: Jackson et al. 2009; 
Lcvrard ct al. 2009; Chang et al. 2010), following 
on earlier work that addressed planet infall even 
during the star's main sequence phase (Gu et al. 
2003; Rasio et al. 1996a, b). Early work on the for- 
mation and migration of giant planets showed that 
migration in the disk could send planets into the 
star (Masset & Papaloizou 2003; Lin & Papaloizou 
1986). The detection of a few unexpectedly close 
massive planets such as WASP-18b (Hellier et al. 
2009) has challenged tidal migration theory to ex- 
plain the statistics of finding planets with such 
small semi-major axes, but it appears that plan- 
ets must continue to migrate towards the star long 
after formation (Hellier et al. 2009), and thus are 
likely to be accreted even after the dissipation of 
the circumstcUar disk. While many other ques- 
tions regarding tidal theory are being addressed, 
it has not been fully considered how much energy 
these massive planets will deliver to their stars. 
Planet consumption has been considered for caus- 
ing contamination of the star by planetary mate- 
rial (Santos et al. 2010; Sandquist et al. 1998; 
Sandquist et al. 2002), even including planet con- 



tamination of white dwarfs (Jura ct al. 2009), and 
the effect of stellar spinup have been considered 
(Pont 2009; Livio & Soker 2002). Planet con- 
sumption during late-stage stellar expansion has 
been widely considered (Soker & Tylenda 2003), 
including as an explanation for the outburst of 
V838 Monocerotis (hereafter V838 Mon; Rctter & 
Marom 2003) and for features of planetary nebulae 
(Soker 1996). These works included consideration 
of the large energy of such consumption (Siess & 
Livio 1999), and work on the energy release has 
been done for the case of accretion onto a white 
dwarf (Shaviv & Starrficld 1988). While this work 
has been done on energy release from the planet 
occurs, it is not clear how planet accretion energy 
will be radiated in most of the many cases. 

Much work has been done on the possibility 
that planets have contaminated stars (Israelian 
2001, 2004; Ashwell et al. 2005) and even white 
dwarfs (Farihi et al. 2010a; Farihi et al. 2010b; 
Jura et al. 2009a; Jura et al. 2009b) with Li 
isotopes and heavier elements. It is important to 
follow through the stages of the end states of plan- 
ets. 

Planet destruction has long been considered to 
occur when a star expands (Shaviv & Starrfield 
1988). It has also been proposed that the bright- 
enings of V838 Mon and a small number of similar 
stars in nearby galaxies might be explained by the 
consumption of planets (Soker & Tylenda 2003). 
The controversy over these events may illustrate 
the difficulty in confirming a particular brighten- 
ing to be a planet destruction event. V838 had 
three brightcnings requiring either the consump- 
tion of three planets, or three separate times of 
mass loss for a single planet. If V383 Mon was 
a planet consumption event, wc; would cixpcct a 
significant fraction of similar brightcnings to be 
explainable by the consumption of single planets 
of different masses. The limited time-series cover- 
age of other candidate events M31-RV and V4332 
Sgr (Bond & Siegel 2006) illustrate the importance 
of having improved data from upcoming time sur- 
veys with foUowup programs. We are not focused 
on planet consumption due to the star's evolution 
due to the previous work on planet consumption 
by evolved stars. However, planet consumption 
that is due as much to the star expanding as to 
tidal migration may be much more likely to avoid 
RLOF and produce a bright event, when the outer 
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layers of the star expand to distances at which the 
planet does not RLOF. 

In Section 2 we address the patterns that lead 
us to consider planet destruction as an important 
event (section 2.1), and given these patterns how 
often it may occur (section 2.2). In Section 3 we 
emphasize the large amount of energy that can 
be released, and address the major question of 
how much of this energy will be released during 
RLOF or collision. In Section 4 we discuss how 
the evolution of tidal migration. In Section 4.1 
we describe the star's rotation as the controlling 
factor allowing stellar tides to cause tidal inspi- 
ral, but in Section 4.2 we consider that magnetic 
braking can slow stellar rotation. We describe how 
RLOF can change the inward migration to alter- 
nating inward/outward migration. In Section 5 
we briefly discuss the inspiral of smaller bodies 
during the planet formation phase. In Section 6 
we advocate future work that can be done dur- 
ing already planned transient astronomy observa- 
tions that may find these events. In Section 7 we 
conclude with an emphasis on observing not only 
the energy release but also seeking to observe the 
material coming off planets that will be a sign of 
planet-destruction related RLOF. 

2. Patterns in planet destruction 

In section 2.1 we address how the distribution 

of planets may be better explained by having lost 
planets to inspiral into the star, and in section 2.2 
we comment on how a luminous signature might 
be rare enough to not have been previously iden- 
tified, but not so rare that current and upcoming 
transient surveys may catch potential events. 

2.1. Planets gone missing 

A major puzzle has been the apparent "pile-up" 
of planets at periods of two or three days, with 
the plausible explanations for bringing planets in- 
ward including Kozai oscillations, secular migra- 
tion, and planet-planet scattering (Mazeli 2008: 
Wu et al. 2007; Faber et al. 2005). So called 
"Hot Jupiters" have long been thought to have 
migrated to their current positions from a more 
distant formation location. Previous efi'orts to 
solve this puzzle sought explanations that have the 
planet absolutely stopping short of falling into the 
star, and indeed, some of the processes such as 



Kozai oscillations do not continue to operate all 
the way to bringing a planet into a star (Fabrycky 
& Tremaine 2007; Wu et al. 2007), though there is 
no reason other than the small cross section that 
planet-planet scattering could not move a planet 
into a collision orbit with the star. Jackson et al. 
(2009) and Levrard et al. (2009) showed that once 
these other process place planets in this pileup, 
that tides caused by the planet on the star may 
be large enough to continue the inward migration 
until the planet is destroyed either by RLOF or 
collision. The rate of tidal migration depends on 
the poorly known amount of stellar tidal energy 
dissipation, which is represented in inverse form 
by the "tidal quality factor" or "Q" value. The 
question of how much energy is dissipated, leading 
to what value of Q*, has been the subject of much 
research which is leading to a consensus that Q* is 
likely a more complicated function than previously 
assumed. Though these tides are still small in the 
pileup, the tides on the star are such a stiff func- 
tion of semi-major axis that they may help define 
the inner edge of the pileup (Jackson et al. 2009; 
Levard et al. 2009; Ibgui & Burrows 2009). De- 
spite the plausible reasons for migration to slow on 
approach to the star in a manner that may explain 
this pile- up (Wu et al. 2007), these reasons do not 
require that the migration stop. We emphasize 
downward migration of the semi-major axis, but 
note that some processes that affect the eccentric- 
ity such as planet-planet interactions do not nec- 
essarily have a barrier preventing eccentricity to 
evolve so close to 1 that the planet could collide 
with the star, releasing slightly more energy than 
inspiralling planets. (Fabrycky & Tremaine [2007] 
show that Kozai oscillations do have a barrier pre- 
venting collision with the star). We emphasize in- 
spiral into the star, which we assume to be much 
more common. 

It had been assumed that after the circumstel- 
lar disk dissipates, planets were considered to pri- 
marily migrate due to tides raised by the star onto 
the planets. Since the star can raise variable tides 
on planets only if when the planets orbit is eccen- 
tric, it was shown that since the tides circularize 
the planets orbit that the planets will indeed be 
preferentially deposited in orbits of two or three 
days. The assumption that planets forever remain 
in these orbits was challenged by Jackson et al. 
(2008a,b,c,d) and Levrard et al. (2009), who no- 
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ticed that at these close-in orbits, the tides raised 
by the planets had been wrongly neglected. Not 
only do planets in three day orbits raise enough 
tides on the star to cause continued slow migra- 
tion in even day-long close-in orbits, but this tidal 
migration is such a stiff function of semi-major 
axis that the inner edge of this pileup may be de- 
termined by planets migrating out of this pileup 
into the star. Jackson et al. (2009) pointed out 
that it may be that the inner age of this pileup 
may be a function of stellar age, as more of the 
"piled up" planets further out migrate into the 
star. Selection effects and tidal spinup may con- 
found the result as originally presented (Jackson, 
personal communication), so improved evidence of 
a stellar age- dependence of minimum planet semi- 
major axis is worth seeking out. Jackson et al. 
(2009) also list several likely observational conse- 
quences of planets likely migrating into the star, 
including searching for rapidly rotating stars and 
stars contaminated by planetary metallicities. 

2.2. Planets gone missing 

If planet destruction of highly massive planets is 
a highly luminous event, it would be rare enough 
to explain why there are no claims of observing 
planet destruction due to inspiral. The only pos- 
sible claims of observing planet destruction events 
are for when the star expands and consumes the 
planet, as in the stars brightening similar to V838 
Mon (e.g. Retter & Marom, 2003). Combining 
the frequencies of "Very Hot Jupiters" and "Hot 
Jupiters" given by Gould et al. (2006) gives one 
in 210 stars that currently hosts a planet that may 
create such an event, but it is unknown how many 
stars have already lost planets. If the current star 
formation rate of the galaxy is 7 stars per year 
(Diehl et al. 2006), then such an event may oc- 
cur once every 30 years per galaxy, and possibly 
more frequently depending on what fraction of the 
original number of planets the current population 
represents. This is a low enough rate to require 
the monitoring of several nearby galaxies by tran- 
sient sky surveys, something that is being started 
with PTF and Pan-STARRS, and will be done 
with greater coverage with LSST (Beatty & Gaudi 
et al. 2008). 

The key question for observing luminous events 
at these rates is whether the events will be fast 
and luminous, or if the gas will be bled off slowly 



enough to not produce dramatic events. The pos- 
sibilities that will lead to luminous events would 
be if enough of a dense part of the planet survives 
and impacts with the star either during main se- 
quence inspiral or when the star expands, and a 
sufhciently rapid runaway RLOF. A slower RLOF 
could be much less luminous, but could allow for 
more opportunities to find close-in planets releas- 
ing sufficiently large amounts of gas that will make 
these objects prime observing targets. 

3. Significant energy release 

The orbital energies of the most massive plan- 
ets is equal to tens of thousands of years of stellar 
luminosity. We plot the orbital energies (potential 
plus kinetic) of the transiting planets closer than 
0.05 AU with full orbital parameters and temper- 
atures available 2010 August 1 in the Extrasolar 
Planet Encyclopedia (exoplanet.eu). We calculate 
the current and "final" orbital energies by adding 
the orbital kinetic energy to the potential energy 
relative to the stellar photosphere (one stellar ra- 
dius). We define the final energy as where what 
we call the start of "planet destruction" to occur. 
We consider "planet destruction" to be occuring 
either at the start of RLOF or if the planet under- 
goes collision with the photosphere. The energies 
are of the order of 10^^ erg as illustrated in Fig- 
ure [TJ This energy is of the order of a small nova. 
We see in Figure [5] that the final orbital energies 
are not much less than the current energies, mean- 
ing most of this energy is released after the start 
of planet destruction. By dividing these energies 
by the yearly luminosity energies of each host star, 
we obtain Figure [3l which shows that the energy 
available is the equivalent of the order of lO'' to 
lO'^ years of stellar output. Expressed this way, 
we again see in Figure |4] that most of this energy 
is release during the process(es) of planet destruc- 
tion. We also can see that if this energy is released 
quickly, then these events have the potential of be- 
ing detected at great distances. In Figure [S] we 
see what the bolometric magnitude would be if 
the energy at RLOF or collision distances were re- 
leased in one day. A real event could release its 
orbital energy moderately more slowly and thus be 
much less bright than this, but still be observable 
in nearby galaxies by upcoming surveys, especially 
by LSST, which will observe to better than mag- 
nitude 18.5 (Beatty & Gaudi 2008). Less massive 
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planets may slowly release their energy through 
RLOF. The wide range of planet accretion events 
will vary, both with planet mass and with when in 
the stages of planet system evolution the destruc- 
tion occurs. 

These large energies are enough to to create 
very luminous events if the energy is released 
quickly. To demonstrate this, we show in Figure [5] 
the brightness that these events could be if the en- 
ergy were all released evenly through one 24-hour 
day. It is neither expected that the energy will be 
released this fast nor that the energy release would 
necessarily be constant, but this arbitrary choice 
is convenient to think of observationally and to 
scale down from. Even if these optimistic bright- 
nesses are scaled down some, these brightnesses 
raise the possibility that current transient surveys 
such as the Palomar Transient Factory (PTF) and 
Pan-STARRs may catch such events in M-31, and 
the LSST could monitor a few more of the near- 
est Milky- Way size galaxies for such events. The 
main challenge may be the rarity of such events 
(discussed below). We advocate work addressing 
the questions of how transient surveys might iden- 
tify observations of planet destruction in both the 
tidal migration and final destruction stages. Dur- 
ing tidal migration, how is the significant energy 
release of massive planets migrating down many 
stellar radii released? How does disk extinction af- 
fect observations of protoplanet accretion? Final 
destruction will occur through different density- 
dependent channels, from RLOF of small planets 
to orbital merging with the stellar photosphere for 
the most massive planets. The process of inward 
migration will be slowed by the planet synchroniz- 
ing the rotation of the star leaving less difference 
between the rotations of the orbit and the star. 
While the speedup of the star's rotation will be 
reduced or reversed by magnetic braking (Barker 
k Ogilvie 2009; Dobbs-Dixon et al. 2004). Mat- 
sumura et al. (2010) find in their simulation of the 
orbits of transiting planets that even though mag- 
netic braking may, in some cases, slow the evolu- 
tion of these orbits, the changes are small enough 
to not greatly affect results such as how planet 
destruction proceeds or how often planet destruc- 
tion might occur. If the final destruction energy 
is released quickly enough, upcoming surveys may 
catch the brightest of these events even in nearby 
galaxies. The challenge will be to identify when 



and how this higher luminosity will produce an 
observable signature. It will be important to iden- 
tify signatures identifiable to upcoming transient 
surveys. We advocate efforts to model and ob- 
serve what we expect will include the most ener- 
getically dramatic events possible in solar systems 
evolution. 

During the final infall, for the most massive 
planets the energy from tides may also become sig- 
nificant relative to the energy output of the star, 
depending on the rate of tidal energy dissipation 
in the star, which is directly related to the stellar 
tidal quality factor Q*. We calculate the energy 
dissipation to be the energy lost to tides on the 
star, which is the negative of the rotational energy 
change of the system (Murray & Dermott, 1999). 
The change of energy of the system is the change 
in the orbital energy minus the energy input into 
speeding up the rotation of the star. We neglect 
the small amount of energy input into speeding up 
the rotation of the planet. 

We relate this energy change to the change in 
semi-major axis. The equations of the change in 
semi-major axis that we use include the change 
of rotation, and are also important in Section 4.1 
where we address how the tidal migration rate is 
coupled to stellar spin-up. To do this we follow 
Levard et al. (2009) in using Neron de Surgy & 
Laskar's (1997) expression that gives the change 
in semimajor axis to second order in eccentricity 
for the evolution of the semi-major axis, assuming 
zero planetary obliquity. 
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where and LOp are the stellar and planetary ro- 
tation velocities respectively, e is the stellar obliq- 
uity, a is the semimajor axis, e is the eccenctricity, 
and n is the orbital mean motion. Q'^ and Q'p are 
the ratios between the present annual stellar and 
planetary tidal quality factors and Qp, and the 
tidal Love number of degree 2 A:2,* and k2^p respec- 
tively. The first term, with Q'^ in the denomina- 
tor, is due to tides within the star, and the second 
term, with Q'p in the denominator, is due to tides 
within the planet. We note that this is similar but 
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not identical to the equation given by Jackson et 
al. (2009). 

The evolution of the semi-major axis is coupled 
to the evolution of the eccentricity, for which we 
use the equation of Jackson et al. (2009), 
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To find the energy input into the star, we follow 
the derivation by Murray & Dermott (1999) that 
this change in energy as a function of the change in 
semi-major axis can be expressed as the change in 
the orbital energy dEorb/dt minus the increase in 
the rotation energy of the star dE,:^rot, and where 
we ignore the much smaller increase in rotational 
energy of the planet as the orbital frequency in- 
creases, 
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The result can be simplified by expressing it in 
terms of the mean motion. 
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where in the case of an inwardly migrating planet 
spinning up a star we have da/dt and [uj — n) less 
than zero. This loss in rotation energy dE/dt goes 
into the star. We show in Figure [B] the amount of 
energy input into the star in units of stellar energy, 
using as an example the XO-3 system in which the 
mass of X0-3b is 11.8 Mj. We show that even for 
a stellar value of 10^, during the last stages of 
tidal migration for the most massive planets the 
energy input into the star could rival the luminos- 
ity of the star. However, if the recent consensus 
(Lanza 2010; Barker & Ogilvie 2010) that for 
F stars such as XO-3 (a type F5V star) is 10^ or 
higher is correct, then the tidal input energy may 



rise to a less significant level. We will discuss the 
question ofthe whether rotational synchronization 
will end migration or if magnetic braking will slow 
the star to continue the planets' migration. We 
leave for further work the important question of 
how fast this energy may be radiated away. 

This paper focuses on actual observations of 
planet destruction due to tidal migration. We 
first briefly comment on previous work on evi- 
dence for planet destruction having already taken 
place, and on previous work on stars expanding 
and consuming their planets as a standard part 
of stellar evolution. This paper empahsizes mas- 
sive planets tidally infalling after planet formation 
is complete, but we end with some comments on 
planet destruction while the circumstellar disk is 
still present. 

4. Tidal migration before and during de- 
struction 

4.1. Tidal migration rate coupled to stel- 
lar spin-up 

Migration into the star is coupled to the ro- 
tational spin-up of the star by the planet, so we 
need to follow and include the rotation of the star, 
which is why we use the form for da/dt (Equa- 
tion [J) of Neron de Surgy & Laskar's (1997). 

The slowdown of the rotation of the star as it 
sheds angular momenta from the system through 
"magnetic braking" is an important question for 
further work. Equation ( [1]) shows that tidal mi- 
gration is much more dependent on the angular 
motion of the planet relative to the star than on 
the planet's absolute rotational motion. The most 
massive planets have sufficient angular momen- 
tum to synchronize the rotation of the star, and in 
some cases (e.g. r Boo) apparently may have done 
so. However, magnetic braking is expected to play 
such a significant role that even all synchronized 
star-planet systems may be expected to be unsta- 
ble to merging, though how fast this will occur 
is of current theoretical work (Barker & Ogilvie 
2007). 

Inward migration will be paused by RLOF, 
which causes outward migration (Chang et al. 
2010; Gu et al. 2003; Paczynski 1966). Though 
estimates have been made (e.g. Gu et al. 2003), 
the equation of state of how much the inner planet 
will swell due to removal of mass from the outside 
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remains a less fully studied but important sub- 
ject for planetary evolution. A delayed swelling 
of the inner planet could lead to sustained RLOF 
even after the planet has started to migrate out- 
ward, leading to further outward migration. It 
will be important to look for signs alternating out- 
ward/inward patterns of migration (Chang et al. 
2010) such as RLOF gas. It is important to de- 
termine how far the outward migration continues 
in order to determine whether currently known 
close-in planets may have already started down 
the planet destruction path. How strongly mag- 
netic braking slows rotation is further complica- 
tion by how the presence of the planet can disrupt 
the steady stellar wind, thus reducing the rate of 
angular momentum loss due to magnetic braking 
(Lanza 2010). 

4.2. Will migration for supermassive plan- 
ets be stopped by synchronization? 

Levrard et al. (2009) raised the question of how 
many planets are unstable to destruction versus 
how many planets will first synchronize the ro- 
tation of the star. Pont (2009) and others have 
begun the study of whether a correlation between 
planet parameters and stellar rotation that might 
indicate such a spin up. In part because of uncer- 
tainties in determining stellar ages the consensus is 
that the question be studied further. Previous age 
determinations have used the stellar slow down in 
rotation, but this cannot be used for planet host- 
ing stars because the rotation rate may have been 
altered by the planet. 

To illustrate the evolution of a very high-mass 
planet, we use X0-3b as a good example system 
because its high mass of 11.8 Mj (Winn et al. 
2008) can be used to illustrate the boundary be- 
tween whether a massive planet is or is not stable 
to infall. Figure[7]shows what would happen to the 
X0-3b system (with period of 3.3 days [Winn et al. 
2008]) following Equation[T]if the star were to now 
be started with zero rotation, where the frequency 
of the planet relative to the star is shown as a 
solid line because it is what determines the migra- 
tion. This frequency is the difference between the 
orbital frequency (dashed line) and w^cose, the 
component the star's rotational frequency (dotted 
line) in the plane of the orbit. This illustrates that 
for such a massive planet that either the star's 
rotation would have to be very small or the or- 



bit would have to be retrograde to not predict 
synchronization based on angular momenta alone. 
Rossiter-McLaughlin measurements have in fact 
found planets in retrograde orbits that would sat- 
isfy this condition. The recent finding that many 
planets, especially the most massive planet sys- 
tems likely including the XO-3 system (Winn et 
al. 2008), have retrograde orbits means that there 
may be more massive planets that collide into their 
stars rather become rotationally synchronized. It 
is clear from Figure [7] that the contribution of the 
inclination is likely the determining factor deter- 
mining the eventual fate of these planets. How- 
ever, even if XO-3 were to have a prograde orbit, 
we can see that the region of stability afforded by 
synchronization could be ended if magnetic brak- 
ing were to slow the system's rotation to the point 
that the period goes under roughly 2 days. At 
this point it reaches a region in which the planet 
is so close to the star that its angular momentum 
is too low to significantly spin up the star. This 
coincides with the stellar tides becoming rapidly 
stronger, as shown by Figure HI which shows the 
time evolution of the semimajor axis of the XO-3 
system in which we have used Equation 1 for zero 
initial stellar rotation. We see a final rapid inspiral 
due to the evolution equations being such a stiff 
function of semimajor axis. Any planet that gets 
this close to a star will rapidly spiral into the star, 
producing a luminous event. 

4.3. Stellar rotation slowed by magnetic 
braking. 

Rotational synchronization becomes unstable if 
the planet's semi- major axis becomes too small, 
such as due to magnetic braking, as illustrated in 
the quick run-up of rotational frequencies in Fig- 
ure [7l Figure [71 which shows the evolution of the 
semi-major axis as a function of time for the same 
parameters as Figure [71 We see that at too small 
of stellar radii, roughly below three stellar radii in 
this case, the planet simply has too little angular 
momenta to continue to spin up the star. A planet 
that has survived to this close distance is likely to 
see its orbit rapidly deteriorate through tidal mi- 
gration, quickly "falling" into the star. Hence, it 
may be unlikely to find planets this close by usual 
planet-search methods due to the short amount of 
time planet may spend there, but the rapid energy 
release may provide a means of finding planets in a 
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final doomed end state. Being so deep in the star's 
potential well (with consequent high orbital veloc- 
ity) such a planet will release a very large fraction 
of the very large energies previously mentioned. 

4.4. RLOF and migration. 

The RLOF transfer of mass from the planet 

to the star causes outward migration of the 
planet due to conservation of angular momentum. 
Paczynski (1966) point out that the separation be- 
tween two orbiting masses is smallest when their 
masses are equal. We follow Paczynski (1966) to 
show that the conservation of angular momentum 
requires that mass transfer from the planet to the 
star results in a larger semi-major axis. We show 
the total angular momentum squared, J^, of the 
planet-star system can be obtained relative to ei- 
ther body by the product of the velocity squared 
of the other body 



Vo = 



GMf 



a{Mi + M2) 



(6) 



multiplied by the mass squared of the first body, 
Ml multiplied by the separation, or semimajor 
axis, squared, a^. Regardless of how we assign 
and Mp to Mi or M2, we obtain. 



=Ga 



[MMpf 

M* + Mr, 



(7) 



where G is the constant of gravitation. When an- 
gular momentum and mass are conserved in the 
system, we have 



a {M,Mp) = 



^-''^^^±^= const. (8) 



We solve for a small change in semimajor axis 
6a resulting from a small transfer of positive mass 

gain to the star of SM^, = 5m from a mass loss 
from the planet of 5Mp = —dm. We obtain. 



6a = 



2a(M* - Mp)6m 



(9) 



M^Mp 

which in the limit where M^ + Mp w M* becomes, 
2a 6m 



6a 



Mp ' 



(10) 



hence RLOF causes outward migration. 

A major unanswered question is whether plan- 
ets will continue to migrate inward as they RLOF 
until they are completely destroyed, or if the an- 
gular momentum transferred through RLOF will 
stall the inward migration and move them back 
outward, and if so, for how long? Will the RLOF 
merely pause the migration, or will a delayed re- 
bound of the planet to the loss of mass move the 
planet outward? How many such inward/outward 
cycles occur before destruction? How much more 
rapid is the outward migration than the inward 
imgration? Might the planets known to be closest 
to their stars have already undergone an RLOF 
and outward migration event, or might we find 
one soon? In other words, what fraction of plan- 
ets are undergoing this event? Gu et al. (2003) 
developed important theory that will be needed 
to address these questions. A very massive planet 
may experience substantial changes in its equation 
of state through either a prolonged pause or pro- 
longed cycles of alternating inward/outward mi- 
gration (Chang et al. 2010). We suggest that a de- 
layed RLOF could maintain an orbital eccentricity 
if more RLOF occurs after periastron than occurs 
before periastron, with the released gas giving the 
planet angular momentum transfer that is unbal- 
anced through its orbit. Such a maintained eccen- 
tricity could enhance tidal heating of the planet, 
further increasing runaway RLOF. Perhaps run- 
away RLOF is sometimes a relatively fast event 
producing short-lived disks of expelled gas, if the 
runaway is not dampened by rapid outward migra- 
tion. Runaway RLOF may be a very rapid event 
(Chang et al. 2010) leading to rapid energy re- 
lease. 

The structure of the released gas may be in- 
clude a disk (or ring) around the star, and/or 
a funnel from the planet towards the star de- 
pending on whether the mass release is enough 
to overcome the stellar magnetic fields (Chang et 
al. 2010). Planet-produced disks of gas should 
be searched for structure, time dependent change, 
and asymmetry. An asymmetric disk (or ring) 
could be a sign of structure that could contribute 
to more massive planets having higher eccentric- 
ities. Rocky planets may respond differently, 
given that the RLOF distance decreases with in- 
creasing density. Similarly, when a massive gas 
planet's higher density core is exposed, the re- 
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maining planet may migrate further in. 

A particularly unstable situation may occur for 
planets of several Jupiter masses, in the regime 
where an increase in mass actually reduces the 
planet's radius. Removal of outer material may 
lead to runaway RLOF as each layer of material 
removed reduces the RLOF radius of the planet 
leading to even more RLOF, until the planet has 
either migrated much farther out or the planet's 
mass is reduced to where more mass loss reduces 
its radius. Gu et al. (2003) describe how adiabatic 
expansion continues to send material out mostly 
through the inner Lagrange point LI. Gu et al. 
(2003) find that RLOF through the outer La- 
grange point L2 is not necessarily exluded, though 
Li et al. (2009) find in their model of WASP- 
12 sufficient loss of material throng LI such that 
material is not lost through L2. Such RLOF can 
cause rapid outward migration of the planet. We 
have shown that donservation of angular momen- 
tum dictates that loss of mass towards the star will 
lengthen the semi-major axis (Paczynski 1966). 
This outward migration brings the Roche distance 
further from the planet to the point that without 
sufficient explansion of the planet the RLOF and 
outward migration will cease, allowing the (slower) 
inward tidal migration to resume. (If a circum- 
stellar disk has not yet fully dissipated, the planet 
could be migrated out into the disk, or the disk 
could be otherwise affected.) Such an effect could 
even contribute to the pileup in planet semi-major 
axes, if some of the Jupiter-mass planets in the 
pileup are actually remnants of planets of many 
Jupiter masses that have already migrated in and 
out one or more times. The star may be signifi- 
cantly brightened by the gas released during such 
an event, possibly causing some of the gas to be ex- 
pelled outward and carry angular momentum out 
of the system. Such an event may have observable 
signatures, either through luminosity variations or 
by signatures of a release of gas. 

4.5. Light curve signature of planet de- 
struction 

Because planet destructio n has only b een mod- 



Siess fc Livio 



eled f or limited case s (e.g., 

" l999ala[ 



Shaviv fc Starrfield 



Soker fc Tvlenda 



1988 : 

20031 ). further work needs to be undertaken to 
identify the light curve that would be produced 
from stellar brightenings due to planet destruc- 



tion. The results of these models shows that cases 
of the planet colliding into the star that the rel- 
evant time scales are on the order of one day, 
leading to our supposing what would be the ob- 
servability of a planet destruction event if the or- 
bital energy would be released in 24 hours. After 
considering the steps that planet destructon will 
go through, we conclude that the light curves may 
change dramatically with planet mass and other 
parameters. The light curve will follow from the 
steps of planet destruction, which may be summa- 
rized as follows: 

• Input of tidal energy into star. 

• RLOF material falling into star. 

• Collision with photosphere will occur for 
more dense planets or cores when the Roche 
distance is less than the stellar photosphere. 

• Movement of material within star, which 
may vary from settling of material in the 
star, to continued orbit within the star (as in 
the case of expanded stars modeled by Siess 
& Livio [1999a]). This will release additional 
energy not included in calculations here. 

Whether and how each step proceeds is highly 
mass dependent, and will also depend on the age 
and type of the star. Other factors including ab- 
sorption by RLOF or whether RLOF affects the 
eccentricity need to be modeled. 

5. Infall of planetesimals during formation 

Planet formation models start with a higher 
number of planetesimals than planets in com- 
pleted systems, and these planetesimals are pre- 
sumed to frequently fall into the star (Trilling 
1998; Chambers 2006). The energy available from 
planetesimal destruction by infall into the star is 
not nearly as great for much smaller objects, but 
the events are likely to be much more frequent, 
and recording such events could provide important 
constraints on planet formation models. Again, 
the important question is whether the energy re- 
lease of such events might be diluted by slow 
RLOF, but perhaps the larger protostar might 
frictionally degrade the orbit to cause a more rapid 
infall to deeper into the protostar. RLOF is an 
important event, producing new disks, funnels, or 
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perhaps strengthening the formation of jets of ma- 
terial from the protostar. The main question may 
be whether the presence of a disk might conceal 
either the energy or mass release cither directly, 
or by how the fluctuations quantity of disk mate- 
rial changes both the environment near the star 
and the brightnesses of young stars. A large frac- 
tion of the photometric observations of accretion 
onto young stars by CoRoT (Alencar et al. 2010) 
show a regularity in their fluctuations that it is 
easy to imagine could be disrupted by a falling 
protoplanet. In looking for exoplanetary system 
obscrvables, it is important to learn whether plan- 
etesimal destruction produces observational signa- 
tures distinct from usual variations in the flow of 
disk material onto the star. Outburst behavior 
from the disk has long been studied (e.g., Zhu et 
al. 2010). Differentiating small protoplanet in- 
fall from such outbursts may be challenging since 
protoplanetary disk systems can have mass infall 
rates of 10~^ to 1O~^M0 per year. Table 1 shows 
how much energy could be released using typical 
solar system objects, for an assumed young sun 
radius of two solar radii. 

6. Future Work 

The full processes of planet destruction as a 
function of the properties of the planet and star re- 
main to be modeled. It can be cixpectcd that the 
mass and composition of planets, and the mass 
and age of the star could all lead to different final 
outcomes for planets. It is important to model 
what happens to the migration of the orbit, how 
the planet's material is removed or if there is a 
collision, and what the effects are on the star. 
Though migration of the semi-major axis, along 
with consumption by evolving stars, will be the 
major routes to destroy planets, it can be asked 
if there will be planets made to collide with their 
stars by evolution in eccentricity or from planet- 
planet scattering. The presence of a disk, and pos- 
sibly other planets may affect this final and most 
active stage of a planet's evolution. Other plan- 
ets and any disk may be affected by the energy 
output. 

Systems undergoing RLOF should be moni- 
tored for variation in luminosity, because RLOF 
processes may likely not be smooth. The gas will 
likely produce periodic negative and positive con- 



tributions to the system's luminosity. With plan- 
ets possessing tens of thousands of years of stel- 
lar luminosity, a sufficiently rapid runaway mass 
loss could rival the star's luminosity. In the case 
of the WASP-12 system, Li et al. (2010) calcu- 
late that gas from the planet may intercept 0.1 
of the stellar visual luminosity, based on work on 
T Tauri disks (Adams et al. 1998; Hartmann et 
al. 1998). Based on this T Tauri disk work, Li et 
al. also report that viscous stress generates heat 
at a rate of ~ 10~^Lo. We propose that observa- 
tions be made not only to detect periodic negative 
or positive contributions to luminosity, but also to 
seek variations in what may be a turbulent pro- 
cess. This is yet one more reason to monitor tran- 
siting planet systems for variations in luminosity. 
We propose full-orbit continuous observations that 
at minimum could elucidate the structure of the 
released material and could provide atmosphere 
information on such misshapen planets, but may 
also reveal whether the release is a turbulent pro- 
cess. Repeated observations near transit, already 
worthwhile to seek transit timing variations, may 
also show variation in transit depth. 

While this paper has been in preparation, much 
progress has been made in modeling (Gu et al. 
2003; Trilling et al. 1998; Li et al. 2010; Lai et 
al. 2010) and observing (Fossati et al. 2010) the 
RLOF of WASP-12b. The models showed that the 
gas detected away from the planet is not merely 
being radiatively removed as is the case for HD 
209458b, but is actually coming off through RLOF 
through an LI Lagrange point nozzle. Chang et 
al. (2010) have modeled a giant planet under- 
going RLOF in the presence of a protoplanetary 
disk. The obvious next step is to extend these 
models to planets similar to close-in planets being 
found by current transit surveys that do not have 
known protoplanetary disks. Understanding the 
response of planets' internal equation of state to 
the removal of outer layers is essential to knowing 
how angular momentum transfer changes the mi- 
gration. Future statistics of the migration status 
of close-in planets will provide a handle on under- 
standing planetary interior equation of states just 
as current work is placing limits on the values of 
planetary and stellar Q* values (Jackson et. al. 
2010, in preparation; Barnes & Taylor 2010, in 
preparation). It will also be important to couple 
successful stastical models of parameter distribu- 
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tion functions (e.g., Jiang ct al. 2010; Jiang ct 
al. 2009; Jiang, et al. 2007) with events such as 
planet destruction that determine these patterns. 
It is to be expected that there will soon be many 
general models describing these events given the 
current high increase in the statistics of such close- 
in planets. 

7. Conclusions 

Searches for transiting planets are now finding 
close-in planets that have probably already begun 
the process of being destroyed. When we observe 
these planets we are seeing tidal migration in its 
most active state, with planets subjected to the; 
highest force and subjected to the fastest sustained 
change. By observing planets in their final stage 
we will learn much about the properties of plan- 
ets. The mass and energy releases accompanying 
planet destruction 

We present the energies available to planets if 
their material is taken from RLOF distance to 
the surface of the star, and demonstrate the sig- 
nificance of this energy by comparing it to the 
amount of time for the star to radiate an equiva- 
lent amount of energy. We note that planet de- 
struction energies compare with the energies of 
nova, and that nova are observed even from nearby 
galaxies. We promote how planet destruction 
should be considered a potential object of inter- 
est for upcoming transient surveys. Because these 
surveys already plan to do the ideal sort of obser- 
vations to find large energy releases, it is impor- 
tant to characterize how this energy release may 
occur. We stress how that this energy release may 
occur at an unknown rate. RLOF may occur fast 
or slow, and there are cases where a planet may 
come in contact with the stellar photosphere, giv- 
ing the fast energy release of a "collision" . 

Planet destruction, however it occurs, offers a 
new phase space of opportunity for the study of 
planets. The manner of the energy release, as 
well as the release of planetary matter, will allow 
planets and stars to be observed in a manner that 
will illustrate properties that we could study in no 
other way. 
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Fig. 1. — Log of orbital energy (kinetic plus po- 
tential) from current orbital semimajor axis of the 
transiting planets closer than 0.05 AU with com- 
plete parameters (listed in the Exoplanet Encyclo- 
pedia [exoplanet.eu] as of 2010 August). 
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Fig. 2. — Log of orbit energy from point of planet 
destruction. This is defined to be when RLOF 
stars or when planet collides with stellar photo- 
sphere. Comparison with Fig refergsstart shows 
that most of the energy is released after the planet 
has started to be destroyed. 
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Fig. 3. — Log of "equivalent years of stellar en- 
ergy" from current orbit. This is the number of 
years of each star's luminosity that is equivalent 
to the orbital energy. 
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Fig. 4. — Log of equivalent years of stellar energy 
from point of destruction as defined in caption to 
Figure [2] Comparison with Figure |3] shows that 
again, after the energies have been rescaled ac- 
cording to stellar luminosities, that most of the 
energy is released after the planet has started to 
be destroyed. 
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Fig. 5. — The bolometric luminosity for an exam- 
ple one 24-hour day release of the orbital energy. 
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Fig. 6. — Energy output of the tides on the star 
during evolution of the semi-major for X0-3b for 
= 10'' and = 10^. Note the time will scales 
linearly with between values of Q*. 
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Fig. 7. — As the planet spirals in, angular momen- 
tum from the star is transferred from the planet 
to the star. We show what the evolution of the 
XO-3 system if the stellar rotation were initially 
at zero, by displaying the rotational frequency of 
the planet with a dashed line and the rotational 
frequency of the star as a dotted line. The all- 
important difference between the two, which can 
be thought of as the planet's rotational frequency 
in the star's reference frame, is shown as the solid 
line. If the star's initial rotation is high enough 
to make this line go to zero, as might be the case, 
the planet is then synchronized with the star. 
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Fig. 8. — Evolution of the semi-major axis given 

in stellar radii for X0-3b for Q, = 10^; the time 
will scale linearly with other values of Q*. The 
planet will inspiral much faster at 4 stellar radii 
and less, a distance at which the planet's angular 
momentum becomes too low to spin up the star 
enough to stop the inspiral. 
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Table 1: Energy release at RLOF distance rela- 
tive to 2 solar radii for representative solar system 
bodies. The change in magnitude relative to the 
sun for a one-day release of this energy is listed, 
followed by the absolute magnitude that such a 
24-hour release of energy would give the sun. 



Object 


Mass 






Radius 


Energy 


Energy 


Change in 


Absolute 




Mq 








Rq 


Years of 


Days of 


Magnitude 


Magnitude 














Lq 


Lq 






Mercury 


1.7 X 10-^ 


3, 


,5 


X 


10-3 


1.9 


480 


6.7 


-1.9 


Venus 


2.5 X 10-^ 


8, 


,7 


X 


10-3 


28 


7000 


9.6 


-4.9 


Earth 


3.0 X 10-^ 


9, 


,2 


X 


10-3 


34 


8600 


9.8 


-5.1 


Mars 


3.3 X lO-'^ 


4, 


,9 


X 


10-3 


3.3 


930 


7.4 


-2.7 


Jupiter 


9.6 X 10-* 


1, 


.0 


X 


10-1 


6600 


2.7 X 10^ 


16.1 


-11.3 


Saturn 


2.9 X 10-4 


8, 


.7 


X 


10-2 


1600 


8.2 X 10^ 


14.8 


-10.0 


Uranus 


4.4 X 10-5 


3, 


.7 


X 


10-2 


300 


1.2 X 10^ 


12.7 


-8.0 


Neptune 


5.2 X 10-5 


3, 


.6 


X 


10-2 


390 


1.5 X 10^ 


12.9 


-8.2 


Pluto 


6.6 X 10-^ 


1, 


.7 


X 


10-3 


0.1 


19 


3.2 


1.5 


Moon 


3.6 X 10-^ 


2, 


.5 


X 


10-3 


0.4 


100 


5.1 


-0.3 


Ganymede 


7.5 X 10-* 


3, 


.8 


X 


10-3 


0.6 


210 


5.8 


-1.1 


Europa 


2.4 X 10-* 


2, 


.3 


X 


10-3 


0.2 


69 


4.6 


0.1 


Titan 


6.9 X 10-* 


3, 


.7 


X 


10-3 


0.5 


190 


5.7 


-1.0 


Rhea 


1.2 X 10-^ 


1, 


.1 


X 


10-3 


8 X 10-3 


3 


1.6 


3.2 


Enceladus 


6.0 X IQ-" 


3, 


.7 


X 


10-4 


4 X 10-4 


0.16 


0.2 


4.6 


Ceres 


4.8 X IQ-io 


7, 


.3 


X 


10-4 


4 X 10-3 


1.4 


0.9 


3.8 



Note. — Table [T] The amount of energy that would be released for gas falling from RLOF to a typical protestellar radius. 
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